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Constitutive Modeling of Crosslinked Nanotube Materials 

G.M. Odegard, S.J.V. Frankland 1 , M.N. Herzog*, T.S. Gates § , C.C. Fay** 

NASA Langley Research Center, Hampton, VA 23681 

A non-linear, continuum-based constitutive model is developed for carbon nanotube 
materials in which bundles of aligned carbon nanotubes have varying amounts of crosslinks 
between the nanotubes. The model accounts for the non-linear elastic constitutive behavior 
of the material in terms of strain, and is developed using a thermodynamic energy approach. 

The model is used to examine the effect of the crosslinking on the overall mechanical 
properties of variations of the crosslinked carbon nanotube material with varying degrees of 
crosslinking. It is shown that the presence of the crosslinks has significant effects on the 
mechanical properties of the carbon nanotube materials. An increase in the transverse 
shear properties is observed when the nanotubes are crosslinked. However, this increase is 
accompanied by a decrease in axial mechanical properties of the nanotube material upon 
crosslinking. 


I. Introduction 

C arbon nanotube (CNT) composite materials have the potential to provide order-of-magnitude increases in 
stiffness-to-weight and strength-to-weight ratios relative to current materials used for aerospace structural 
applications. These properties are especially important in the design and development of ultra-light-weight 
aircraft such as new classes of Unmanned Aerial Vehicles (UAVs), in which the primary requirements are those 
associated with long-duration, high-altitude flights. To facilitate the development of CNT-based materials for this 
purpose, constitutive relationships must be developed that predict the bulk mechanical properties of the materials as 
a function of the molecular structure. 

Within the past few years, considerable effort has been expended to synthesize CNT/polymer composites that 
take advantage of CNT properties by enhancing the load transfer between the CNTs and the adjacent polymer 
molecules. 1 ' 9 One approach to achieving this is to form chemical bonds between CNTs and adjacent polymer 
molecules (functionalization). Despite the potential increase in load transfer that follows from functionalization 
with respect to the load transfer that occurs without the presence of a chemical bond, it has been recently 
demonstrated that functionalization itself can measurably affect the structure and mechanical properties of 
CNTs, 1011 which can, in turn, ultimately alter or perhaps degrade the bulk mechanical properties of CNT/polymer 
composite materials. As an alternative approach to increasing load transfer efficiency between the CNT and 
surrounding polymer, NASA Langley Research Center has recently developed a material in which single-walled 
CNTs are covalently bonded together (crosslinked) with the chemical linking agent l,3-bis(4-aminophenoxy-4’- 
benzoyl)benzene (1,3-BABB). In this material, the short, organic linker units allow for direct CNT-to-CNT load 
transfer, therefore, possibly allowing for improved overall mechanical properties with respect to pure CNT-based 
materials and CNT/polymer materials in which the CNTs are functionalized to the polymer only. Since it is 
expected that the bulk mechanical properties in this new class of materials are expected to be affected by intrinsic 
and sometimes subtle changes in the molecular structure, a multi-scale modeling approach must be developed to 
allow for accurate design of materials and provide the constitutive relationships necessary for macro-scale analysis 
methods. 

The objective of the present paper is to develop non-linear constitutive models for these crosslinked CNT 
materials and to examine the influence of crosslink density on the overall, continuum-level elastic properties. 
Following a brief description of the materials, the equivalent-continuum modeling approach is described in detail. 
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The approach includes the determination of the molecular structures of the materials using molecular dynamics 
(MD) simulations, the development of a non-linear elastic, continuum-based constitutive law using a 
thermodynamic energy technique, and the determination of the equivalent-continuum material parameters via energy 
equivalence. The constitutive relationships developed in this paper enables the non-linear elastic behavior of the 
material to be described in terms of applied strain, which is necessary for the full application of MD techniques in 
the equivalent-continuum modeling approach. Results from these analyses are compared, and the effect of the 
crosslinking in these materials is discussed in terms of engineering properties. 


II. Materials 

The constitutive models developed in this study are for four material systems: a bundle of aligned single-walled 
CNTs with no crosslinks (0% crosslink density), two bundles of aligned single-walled CNTs with finite crosslink 
densities (0.3% and 0.45%), and the amorphous pure crosslinking material. The crosslink density is defined as the 
ratio of the number of functionalized carbon atoms in the CNTs (due to crosslinking) to the total number of carbon 
atoms in the CNTs. All of the CNTs are achiral (10,10) nanotubes with corresponding diameters of 16.9 A. The 
crosslinking material is l,3-bis(4-aminophenoxy-4’-benzoyl)benzene (1,3-BABB) (Fig. 1). The simulated form of 
1,3-BABB is modified to replace the -NH 2 groups with hydrogens. Figure 2 shows an example of how the CNTs 
are crosslinked. Details on the processing and characterization of these materials can be found elsewhere. 12 


O O 



Figure 1. CNT crosslinking material 1,3-BABB 


III. Equivalent-Continuum Modeling 

The nonlinear-elastic properties of the four materials systems were determined using an equivalent-continuum 
modeling method. 11 ' 13 ' 14 This approach consisted of three steps. First, representative volume elements (RVE) of the 
molecular structure and equivalent-continuum models were chosen that accurately described the bulk structure of 
the material. Next, a constitutive law that described the behavior of the equivalent-continuum model was 
determined. Finally, the energies of deformation of the two models were equated under identical sets of boundary 
conditions to determine each of the elastic constants in the constitutive equation. Each of these steps is described in 
detail below. 

A. Representative Volume Element 

The RVEs of the molecular and continuum models where selected using the molecular structure of the four 
materials as a guide. The equilibrium molecular structures of each of the four material systems were determined 
using MD simulations. The RVE geometry selected was a parallelepiped box. The MD simulation box was then 
prepared with 9 aligned single-walled CNTs (except for the pure crosslink material) with the proper amount of 
crosslink material and functionalization sites that corresponded to 0.0%, 0.3%, and 0.45% crosslink densities. The 
CNTs were aligned parallel to the x r axis of the model. Periodic boundary conditions were used in the simulations, 
thus the CNTs were modeled as having an infinite length. The bonded and non-bonded interactions of the atoms 
were described by the AMBER force field, 15 without electrostatic interactions. The aromatic rings of the 1,3-BABB 
were treated as rigid bodies in the simulations. The initial configurations for each of the four systems were 
compressed using molecular dynamics to minimize the configurational energy, as the size of the simulation box was 
gradually decreased. An nPT (constant number of atoms, pressure, and temperature) algorithm was subsequently 
applied at a temperature of 298 K and pressure of 1 atm for up to a few picoseconds to establish the equilibrated, 
undeformed molecular structures. The resulting equilibrium molecular structures are shown in Figs. 3-6. DL-POLY 
was used for these simulations. 16 Subsequently, the RVE of the continuum models were also formed as solid 
parallelepipeds with dimensions identical to those of the molecular models, as shown in Figs. 3-6. 
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Single-walled 
carbon nanotube 



Front view 


Side view 


Figure 2. Covalently bonded single-walled CNTs 


B. Constitutive Equation 

Many materials experience a nonlinear-elastic response when subjected to static loads . 17 18 The term nonlinear- 
elasticity can have two different meanings in the description and analysis of material behavior. In the context of this 
study, the material is nonlinear-elastic in the constitutive sense, not the geometric sense. That is, there is a one-to- 
one, nonlinear relationship between applied stresses and strains. This study does not consider possible higher-order 
forms of the stress and strain tensors that reflect the changes in boundary conditions that can occur for large 
deformations. 

A constitutive model that is typically used for nonlinear elastic materials is the Ramberg-Osgood 
formulation , 19 ' 20 which was originally developed for isotropic materials with the nonlinear tern associated with the 
applied stress. However, because the materials considered in this study are anisotropic, and because it is difficult to 
apply force-based boundary conditions in MD simulations, an anisotropic, strain-based, nonlinear-elastic 
constitutive law needed to be established. 

Perhaps the most fundamental approach to developing continuum-based constitutive equations for materials is 
using thermodynamics-based energy functions . 17 This approach allows for the selection of a set of independent 
variables, such as strain and temperature, that are used to describe overall material behavior through a state function. 
Since it is necessary for a constitutive law to be expressed in terms of strain when interfacing directly with MD 
techniques, this thermodynamics-based approach is utilized here to describe the mechanical behavior of crosslinked- 
CNT materials. 
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Molecular model 


Equivalent-continuum 

model 


a = 49.9 A a = 90.0° 

b = 49.6 A (3 = 90.0° 

c = 53.8 A (j) = 59.8° 

Figure 3. RVE of CNT bundle with 0% crosslinking 



Molecular model Equivalent-continuum 

model 


a = 48.1 A 

a = 90.0° 

b = 64.5 A 

(3 = 90.0° 

c = 53.8 A 

tj) = 60.0° 

Figure 4. RVE of CNT bundle with 0.3% crosslinking 
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Molecular model 



Equivalent-continuum 

model 

a = 90.0° 

(3 = 90.0° 

4 = 60.0° 


a = 57.5 A 
b = 56.5 A 
c = 53.8 A 

Figure 5. RVE of CNT bundle with 0.45% crosslinking 



X \ 

Molecular model 



Equivalent-continuum 

model 

a = 90.0° 

(3 = 90.0° 

(|> = 60.0° 


a = 54.2 A 
6 = 86.5 A 
c = 27.6 A 

Figure 6. RVE of pure 1,3-BABB crosslinking material 
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For a uniform, isothermal, homogeneous deformation of a volume of an elastic material, the corresponding 
uniform stress in the volume is 17 


_ 1 dF 

n " f n:„ 


( 1 ) 


where V is the volume of the RVE, s,y is the strain tensor (ij = 1,2,3) and F is the Helmholtz free energy, which is 
henceforth referred to as the strain energy. At this point, it is assumed that the only dependent variables in the strain 
energy are the applied strains, i.e. F = F(s, ; ). The functional form of the strain energy can be further restricted by 
considering the invariance properties of the material such that the strain energy remains invariant with respect to the 
coordinate transformations expressed by the material symmetry. Inspection of the molecular models indicates that 
the constitutive law needs to represent a material with orthotropic symmetry. Therefore, the strain energy must be 
expressed in terms of the invariants of an orthotropic material, which are 21 


2 2 2 

^1P ^22’ ^33’ ^23’ ^13’ ^12’ ^23^13^12 


( 2 ) 


For convenience, and because multiples of invariants are also invariants, the strain energy is expressed as 

—F — f (Ejj, s 22 , s 33 , s 23 , s 13 , e 12 , 6 22 E 33 , S u S 33 , S n S 22 ) 

The functional form of the strain energy is further specified to be sums of individual invariant functions 


where 


y F - f\ ( S U ) + fl ( S 22 ) + f ( e 33 ) + f* ( S 23 ) + fs ( S 13 ) + fb ( E 12 ) 


' fl ( E 22 E 33) + fs ( E 1 1 S 33 ) + f<> ( E 1 1 E 22 ) 


f ( E n ) = ~~rAi | E n | 
v ’ n, + 1 11 


fl ( E 22 ) = ~ 7 ^22 

v ’ n 2 + 1 

A ( 4 ) 

V ’ /?, + 1 




(n,+l) 


A 2 


2 |t("3+1) 


/ 4 ( E 2 3 ) = 


(« 4 +i) 

2" 5 


A \p z 
^44 & 23 


/ 5 (4)=^— ty4 5 | e 

\ n 5 +1 ) 


f 6 ( E 12 ) = 


(« 6 +1 ) 


Ae 4 12 


(n 4 +i) 

2(«6+i) 


fi (s 22 s 33 ) — A 3 s 2 2 s 

fs ( E 1 1 S 33 ) = ^13 E 11 E : 

f 9 ( E 1 1 S 22 ) = ^12 E 11 E : 


33 

33 

22 


( 3 ) 


( 4 ) 


( 5 ) 


( 6 ) 
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In Eqns. (5) - (7), n k (k = 1-6) and Ay are material constants. The power-law relationships in Eqns. (5) and (6) were 
chosen in order to yield a form of the constitutive relation that would accurately represent the material deformation 
yet be simple enough to lend itself to the numerical simulations. The constitutive equation that results from Eqns. 
(1) and (4) - (7) are 


CT 11 — A l^ll "*■ ^13^33 

CJ 22 — A n E n + A 22 S 22 + yT 3 £ 33 

CT 33 — ^ 13^11 ^ 23^22 ^ 33^33 /q-, 

„ (o) 

CT 23 — ^44 7 23 4 

CT 13=4 5 Yl3" 5 
CT 12 — ^66 7 12 

where jy = 2sy. From Eqn. (8) two immediate observations can be made. First, the constants Ay are similar to the 
elastic stiffness tensor components in linear-elastic materials. Second, all six components of stress exhibit some 
element of the power-law response with respect to strain. Therefore, it is noted that for a strictly linear elastic 
material, n k = 1, and for a nonlinear-elastic material, n k < 1. As the material is unloaded, both the strains and stresses 
will return along the original loading path. Upon full unloading, there will be zero residual deformation. 

C. Energy Equivalence 

To determine each of the material parameters in Eq. (8), the energies of deformation of the continuum and 
molecular models were equated for 9 sets of boundary conditions. Each set of boundary conditions corresponded to 
homogeneous deformations in the continuum models. The displacements applied at the boundaries of the RVE are 
generalized by 


u, (. B ) = ZjXj (9) 

where B is the bounding surface of the RVE, Xj is defined in Figs. 3-6, and the summation convention associated 
with repeated indices is used. The 9 sets of boundary conditions and the material parameters determined for each set 
are listed in Tables 1-3. The deformations described in Tables 1-3 correspond to axial, plane-strain bulk, and shear 
deformations, respectively. Unspecified strain components in Tables 1-3 are zero valued. 


Table 1. Boundary conditions for axial deformations 


Materials parameters 

Boundary conditions 

Boundary displacements 

An, n\ 

sn = e 

mi (B) = ex\ 
112(B) = 0 
u:(B) = 0 

An, 77-2 

S 22 = e 

u\(B) = 0 
112(B) = e %2 
m 3 (B) = 0 

A 33 , 773 

II 

CO 

CO 

CO 

ui(B) = 0 
112(B) = 0 
112(B) = exs 
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Table 2. Boundary conditions for plane-strain bulk deformations 


Materials parameters 

Boundary conditions 

Boundary displacements 

A23 

£22 = £33 = e 

u\(B) = 0 
112(B) = exi 
113(B) = ex 3 

Al 3 

£11 “ £33 = e 

u\(B) = ex\ 
112(B) = 0 
113(B) = ex 3 

An 

8n = 822 = e 

u\(B) = ex\ 
112(B) = ex2 
113(B) = 0 


Table 3. Boundary conditions for shear deformations 


Materials parameters 

Boundary conditions 

Boundary displacements 

A 44, fU 

S23 = y/2 

u\(B) = 0 
112(B) = (y/2)x3 
113(B) = (y/2)x 2 

A 55, II5 

S13 = y/2 

u\(B) = (y/2 )X3 
112(B) = 0 
113(B) = (y/2)xi 

^ 66 ? ^6 

so = y/2 

u\(B) = (y/2)x2 
u 2 (B) = (y/2)x\ 
113(B) = 0 


While the strain energies of the continuum models were determined for each set of boundary conditions using 
Eqns. (4) - (7), the energies of deformation of the molecular models were determined using MD simulations. For 
the boundary conditions listed in Tables 1-3, four magnitudes of axial and plane-strain bulk strains were applied; e = 
0.5, 1.0, 1.5, and 2.0%; and four magnitudes of shear strains were applied; y = 0.5, 1.0, 1.5, and 2.0%. For the 
molecular models, each of the deformations were applied to the equilibrium molecular structures shown in Figs. 3-6 
by deforming the MD simulation boxes and all of the atoms in the models according the applied strain field. An 
nVT (constant number of atoms, volume, and temperature) simulation was subsequently used for each deformation 
to allow the box dimensions to remain fixed while the atoms were allowed to move into new equilibrium positions. 
The simulations were run up to 2 ps at 298 K using the TINKER modeling package 2 ~ and the AMBER force field 15 
without electrostatic interactions. The energy of deformations of the molecular models were averaged over the final 
400 fs of each simulation, and the standard error of the energy fluctuations were recorded. 

The force field parameters used in the MD simulations are listed in Tables 4-8. In Tables 4-8, R mi „ and E min are 
the Van der Waals radius and well depth, respectively; K, and r are the bond-stretching constant and equilibrium 
distance, respectively; K 0 and 0 are the bond-angle constant and equilibrium angle, respectively; and N pat h, VJ2, y, 
and n are the number of bond paths, torsion magnitude, phase offset, and periodicity of the torsion, respectively. 
Further details of these parameters may be found elsewhere. 15 The carbon atoms in both the CNT and aromatic 
groups of the crosslinks were modeled as aromatic carbons with sp~ hybridization, denoted by CA in Table 4 (which 
is consistent with the AMBER notation). The carbon atoms in the CNTs that were located where the 
functionalization points exist were modeled as having sp 3 hybridization, denoted as CT in Table 4. The oxygen 
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atoms in the ether groups of the crosslinks were modeled as having sp hybridization and denoted by OS. The 
oxygen and carbon atoms in the carboxyl groups of the crosslinks were modeled as having sp~ hybridization and are 
denoted by O and C, respectively. The hydrogen atoms on the aromatic groups in the crosslinks are denoted by HA. 
Fig. 2 shows an example of how the CNTs were crosslinked in the simulations for this phase of the analysis. 

With the molecular and equivalent-continuum energies in hand, the next step was to determine the material 
parameters of Eq. (8). For each set of axial and shear boundary conditions, the computed energies (per unit volume) 
from the molecular models for all four strain levels were fitted with curves using least squares regression, 23 as a 
function of the square of the applied strain. The material parameters Hu, ri\, A 2 2 , n 2 , A 33 , n 3 , A 44 , n 4 , A 55 , n 5 , A 66 , and 
n h were determined for each material using the curve fits and Eqs. (4) - (6). For the plane-strain bulk boundary 
conditions, the material parameters A 23 , A l3 , and A u were determined using the axial parameters A n , n h A 22 , n 2 , A 33 , 
« 3 and Eqs. (4), (5), and (7). Parameters A 23 , A I3 , and A n were adjusted to minimize the correlation coefficient 23 
between the energies of the molecular and equivalent-continuum models over the entire applied strain range. 


Table 4. Van der Waals parameters 


Atom 

Rmin C^-) 

E min (kcal/mole) 

C 

1.9080 

0.0860 

CA 

1.9080 

0.0860 

CT 

1.9080 

0.1094 

HA 

1.4590 

0.0150 

0 

1.6612 

0.2100 

OS 

1.6837 

0.1700 


Table 5. Bond-stretching parameters 


Bond type 

K t (kcal/mole/ A 2 ) 

r (A) 

C-CA 

469.0 

1.409 

C-0 

570.0 

1.229 

CA-CA 

469.0 

1.400 

CA-CT 

317.0 

1.510 

CA-HA 

367.0 

1.080 

CA-OS 

320.0 

1.410 
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Table 6. Bond-angle parameters 


Bond type 

Kq (kcal/mole/radian") 

0 (degrees) 

CA-C-CA 

63.00 

120.00 

CA-C-0 

80.00 

120.40 

C-CA-CA 

63.00 

120.00 

CA-CA-CA 

63.00 

120.00 

CA-CA-CT 

70.00 

120.00 

CA-CA-HA 

35.00 

120.00 

CA-CA-OS 

50.00 

109.50 

CA-CT-CA 

63.00 

114.00 

CA-OS-CA 

60.00 

109.50 


Table 7. Torsional parameters 


Bond type 

Npath 

V n /2 (kcal/mole) 

y (degrees) 

n 

CA-C-CA-CA 

4 

14.50 

180.0 


CA-C-CA-0 

4 

14.50 

180.0 


C-CA-CA-CA 

4 

14.50 

180.0 


C-CA-CA-HA 

4 

14.50 

180.0 


CA-CA-CA-CA 

4 

14.50 

180.0 


CA-CA-CA-CT 

4 

14.50 

180.0 


CA-CA-CA-HA 

4 

14.50 

180.0 


CA-CA-CA-OS 

4 

14.50 

180.0 


CT-CA-CA-CT 

4 

14.50 

180.0 


CT-CA-CA-HA 

4 

14.50 

180.0 


HA-CA-CA-HA 

4 

14.50 

180.0 


HA-CA-CA-OS 

4 

14.50 

180.0 


CA-CA-CT-CA 

6 

0.00 

0.0 


CA-CA-OS-CA 

3 

1.15 

0.0 
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Table 8. Improper torsional parameters 


Bond type 

Npath 

VJ2 (kcal/mole) 

y (degrees) 

n 

CA-CA-CA-CT 

1 

1.1 

180.0 

2 

CA-CA-CA-HA 

1 

1.1 

180.0 

2 


IV. Results 

The calculated material parameters for all four materials are shown in Table 9. For the CNT bundle systems with 
and without crosslinking, all of the exponent («,) values were 1.0 for the given range of strains (up to 2 %) that were 
applied in the modeling. Thus, for strains under 2%, these systems exhibit linear-elastic behavior. For the same 
strain range, the 1,3 -BABB showed significant nonlinearity under axial deformation, as indicated by the lower 
values of n u n 2 , n 2 . The 1,3-BABB showed moderate nonlinearity when subjected to shear deformation in the x 2 -x 2 
plane (/? 4 ). 


Table 9. Material parameters (A,j in GPa, /?, is unitless) 


Parameter 

0% crosslink 
density 

0.3% crosslink 
density 

0.45% crosslink 
density 

1,3-BABB 

A 11, n\ 

570.0, 1.0 

440.4, 1.0 

463.6, 1.0 

0.3, 0.1 

A22, n 2 

39.5, 1.0 

14.3, 1.0 

4.2, 1.0 

2.7, 0.7 

A33, tl 3 

40.1, 1.0 

3.4, 1.0 

0.8, 1.0 

0.4, 0.1 

A23 

23.7 

4.1 

2.1 

0.1 

A 13 

9.4 

2.8 

6.5 

0.1 

An 

12.1 

1.8 

5.4 

9.0 

A44, ri4 

1.0, 1.0 

28.8, 1.0 

22.1, 1.0 

12.0,0.9 

A 55 , 115 

32.1, 1.0 

27.7, 1.0 

36.5, 1.0 

5.3, 1.0 

^ 66 ? ^6 

38.8, 1.0 

36.4, 1.0 

33.8, 1.0 

9.8, 1.0 


From Table 9, the effects of the crosslinking on the mechanical properties of the CNT bundles can be discerned. 
For A u the stiffness parameter in the nanotube axial direction, a reduction of 23% occurred when the CNTs were 
crosslinked at a crosslink density of 0.3 %. Little change in A n was observed upon further crosslinking to a 
crosslink density of 0.45%. More significant losses were sustained by the transverse stiffness parameters A 2 2 , and 
A 33, with reductions of up to 92% at a crosslink density of 0.3%. Furthermore, A 2 2 and A 22 continued to decrease 
significantly upon further crosslinking to a crosslink density of 0.45%. For A 23 , a total loss of 98% occurred 
between the non-crosslinked and the 0.45% crosslinked systems. Therefore, it is clear that the crosslinking of the 
CNTs degrades the mechanical properties of the CNT bundles when these bundles are deformed parallel to the 
principle axes of the material (Figs. 3-6). 

The material response to plane-strain bulk deformations is evident by noting the calculated values in Table 9. 
The transverse plane-strain bulk deformation parameter, A 2 3 , decreased by 83% and 91% for the 0.3% and 0.45% 


11 

American Institute of Aeronautics and Astronautics 





crosslink densities, respectively. The longitudinal plane-strain bulk parameters, A 13 and A l2 , decreased significantly 
(up to 85%) when the CNTs were crosslinked to a density of 0.3%. However, the same parameters increased (up to 
a factor of 3) upon further crosslinking up to 0.45%. Therefore, while the parameter A 2 3 shows a clear degradation 
upon CNT crosslinking, there was no clear effect of crosslinking on the parameters A\ 3 and A 12 . 

The effects of crosslinking on the shear moduli are also evident from the data in Table 9. Under a transverse 
shear load, the resulting mechanical parameter, A 44 . increased by a factor of 29 when the CNTs were crosslinked up 
to 0.3%. Little change was observed under the same loading condition upon further crosslinking. The response to 
longitudinal deformation (A 55 , A 66 ) was not significantly affected by the crosslinking in the bundle systems. 
Therefore, the crosslinking of the CNT enhanced the transverse shear properties while having little effect on the 
longitudinal shear properties. 

For completeness, the stiffness parameters for the pure crosslinking material 1 ,3-BABB are also included in 
Table 9. As expected for an organic molecular solid, the pure material had relatively low mechanical properties 
compared with the CNT materials. In particular, the axial mechanical parameters A n , A 2 2 , A 33 , n,, n 2 , and n 3 ; the 
longitudinal shear parameters A 55 and A 66 ; and the plane strain bulk parameters A 23 and A 13 were significantly 
smaller for the 1,3-BABB than those of all three of the CNT systems. 


V. Summary 

In this study, a nonlinear constitutive model was developed for CNT-bundle materials in which the CNTs had 
varying amounts of crosslink densities. An equivalent-continuum modeling approach was used to develop the 
continuum-based, non-linear elastic constitutive relationship that incorporated the subtle changes in molecular 
structure upon crosslinking. In the equivalent-continuum modeling, MD simulations were performed to determine 
the equilibrium molecular structure of these systems as well as to determine the molecular structure and 
corresponding energy changes when subjected to applied strains. As a result, the effects of material nonlinearity and 
the influence of the crosslinking on the mechanical properties of the CNT materials were discerned. 

Even though the CNT materials did not exhibit any nonlinear behavior, the pure crosslinking material was 
nonlinear, as indicated by the lower values of the power-law exponents in the constitutive equation. Overall, the 
presence of the crosslinks had significant effects on the mechanical properties of the CNT bundles. In particular, the 
crosslinks degraded the axial mechanical properties and the transverse plane-strain bulk properties of the CNT 
bundles. A decrease of up to 98% and 91% was observed for the axial and transverse plane-strain bulk parameters, 
respectively, when a 0.45% crosslink density was applied to CNT bundles. However, the crosslinking greatly 
enhanced the transverse shear properties. A nearly 30-fold increase was observed for the transverse shear parameter 
between the non-crosslinked and 0.3% crosslink density systems. 
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